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Abstract. We present an overview of the occurrence and properties of atomic gas associated with compact radio 
sources at redshifts up to z — 0.85. Searches for Hi 21cm absorption were made with the Westerbork Synthesis 
Radio Telescope at UHF-high frequencies (725-1200 MHz). Detections were obtained for 19 of the 57 sources 
with usable spectra (33%). We have found a large range in line depths, from r = 0.16 to r < 0.001. There is 
' a substantial variety of line profiles, including Gaussians of less than 10 kms -1 , to more typically 150 kms -1 , 

as well as irregular and multi-peaked absorption profiles, sometimes spanning several hundred kms . Assuming 
uniform coverage of the entire radio source, we obtain column depths of atomic gas between 1 x 10 19 and 3.3 x 10 21 
(T S p/100 K)(l//) cm -2 . There is evidence for significant gas motions, but in contrast to earlier results at low 
^ \ redshift, there are many sources in which the Hi velocity is substantially negative (up to v = —1420 kms -1 ) with 

respect to the optical redshift, suggesting that in these sources the atomic gas, rather than falling into the centre, 
may be be flowing out, interacting with the jets, or rotating around the nucleus. 
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1. Introduction quence of circumnuclear starbursts and accretion onto the 

nuclear black hole (e.g., Sanders et al.[l988). 

The energy output from AGN is commonly agreed to be 

powered by accretion onto a super-massive black hole, re- TT . .. . . . , . . . ., , „ , 

... r . „ „ „ . „ ,. . However, while this picture is becoming well denned 

suiting from the mlall of gas transported from the host . „ P , . . .. ' , .. 

. . , , . . , _ ,, .I,, ,. tor Seyterts, very little is known regarding the situation m 

galaxy down to the central regions (e.g., Rees 1984 . ° ° . , . „, T . 

° , , ,. ,n r ; i — ■ powerful radio galaxies, the mostly radio-quiet AGN m 

Recent detailed and statistical studies of Seyfert galaxies . . . . ,. , , TT ; , , , _ , 

/ ^. , ^ , , irmrm i • • { . .i a large nearby galaxy sample studied by Ho et al. 1997 

(e.g., Cid lernandes et al. 2001 are beginning to chart the ° . , / , , f , . V 



evolutionary sequence that leads to the production of an 



often reside in early- type disk galaxy hosts, many of which 

, , • • , - 1 i , , r , , , have spiral structures, and the amount of gas needed to 

active galactic nuclei m late- type (Sd through SO) galaxies , . , . . ., ,. ., „ , ,. r 

, 1 . „. , , , , . , i , / sustain the low- luminosity activity for the expected hfc- 

and the effect of that evolution on the host galaxy (e.g., .. . ,., ., , , „ , , , . ,. , , .^ AT 

_ . . „ . —-— rT ^. . , n ; • time is readily available. But looking at radio-loud AGN, 

Storchi-Bergmann et al.EDQD. The picture that is emerg- MM ^ ^ m ^ ^ ^ ^ ^ ^ 8Q% rf ^ 



ing is one of a close tie between the galaxy environment 
and galaxy-galaxy interactions leading to a temporal se 



hosts of 46 low-redshift {z < O.f) 3CR radio galaxies are 
ellipticals. At larger redshifts, Dunlop et al. H2003(l have 

shown that the host galaxies of radio-quiet and radio-loud 

Send offprint requests to: R.C. Vermeulen, quasars, as well as radio galaxies, are massive ellipticals. 

rvermeulen@astron . nl These may of course themselves be merger products of late 
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type galaxies. True spiral hosts of powerful radio galaxies 
seem to be very rare (see also Ledlow et al. 1998 ). 

The predominance of elliptical hosts is often thought 
to be puzzling, since these are usually assumed to be gas 
poor. It should be noted first that major galaxy merger 
products are likely places where gas has lost angular mo- 
mentum and can feed into the central engine. Apart from 
that, Walsh et al. (1983 have shown that in ellipticals the 
amount of long-wavelength IR emission, indicative of the 
presence of dust and gas, is correlated with the occurrence 
and strength of a central radio source. More recently, it 
has been shown that early-type galaxies can contain gas 
in both molecular form (see e.g. Rupen [T997; Knapp & 
Rupen 1996), and in atomic form (see e.g. Bregman et al. 
1992 Oosterloo et al. H999|l , and this gas could act as a fuel 
reservoir for the central engine. Indeed, central structures 
of gas are frequently observed in nearby radio sources; for 
instance, 17 of 19 nearby FRI radio galaxies observed with 
the HST by Verdoes Kleijn et al. (1999) have dust struc- 
tures in their centres. Such dust disks are also reported in 
molecular (e.g in 4C12.50, Evans et al. I1999J) and atomic 
lines, e.g. in NGC4261 (Van Langevelde et al. l2HU0|l and 
NGC1052 (Vermeulen et al. l2003al) . Furthermore, obser- 
vations of nearby radio-loud AGN by van Gorkom et al. 
( 198!]}) have directly shown, in about 30% of the galaxies 
studied, the presence of cold gas visible through the Hi 
A = 21 cm line of neutral hydrogen in absorption towards 
the central radio source. In a recent study, Morganti et 
al. H2001f) have found evidence that in their radio galaxy 
sample the H I detections are related to circumnuclear tori. 

At cosmological redshifts, the galactic medium is 
largely unexplored. But the recent upgrade of the 
Westerbork Synthesis Radio Telescopes (WSRT) has 
opened up a new avenue of research. The array of 14 tele- 
scopes, with its new multi-channel correlator (DZB), offers 
good sensitivity as well as a superior capability to reject 
external interference compared to single dish telescopes. 
The wide band UHF-high receivers, spanning 725-1200 
MHz, enable studies of the occurrence and kinematics of 
atomic gas seen by means of the redshifted H I line in ab- 
sorption against radio sources almost out to z = 1. We 
here present the combined results of a number of H I line 
searches performed at the WSRT. These projects had a 
variety of goals, all related to the properties of cold gas 
in compact radio sources with a total extent less than a 
few arcseconds, within the appropriate redshift range, in 
the Northern sky, and with a flux density above a few 
hundred mjy (in order to have a reasonable opacity de- 
tection threshold). A majority of the sources in this pa- 
per are classified as Gigahertz Peaked Spectrum (GPS) or 
Compact Steep Spectrum (CSS) sources, believed to be 
intrinsically compact (e.g. Fanti et al. 1995, O'Dea 1998). 
Further, a subset of the sources in this paper are end-on 
classical doubles (including both quasars and galaxies). 
As a result of coordination between the present authors, 
we here present the complete observational results for the 
sources in the original programmes for which usable spec- 
tra were obtained, plus some statistics on the absorber 



depth and kinematics. In a follow-up paper, using a more 
tightly defined sample of GPS and CSS sources including 
data from the literature, Pihlstrom et al. ( 2003) discuss Hi 
absorption in the context of other properties, and show, 
for example, the existence of an anti-correlation between 
absorption depth and linear size of the radio source. 

2. WSRT observations and data analysis 

This paper presents the combined results of a number of 
different WSRT projects, performed between 1997 and 
2001, using the UHF-high receivers. While there were 
slight variations in the observational and data reduction 
methods used, the description below adequately describes 
the essentials. The data reduction was carried out using 
the NRAO AIPS and Caltech DIFMAP packages, typi- 
cally switching back-and-forth to make use of their respec- 
tive strong points. 

Typically, each target was first observed for a few 
hours, in each of two orthogonal linear polarisations, with 
128 spectral channels covering a 10 MHz wide observing 
band (the maximum available when the surveys started), 
centred at the frequencies predicted for Hi based on the 
optical redshifts (these have some uncertainty, as dis- 
cussed in Section|2J). Doppler tracking was used; it ensures 
that spectra taken separately at arbitrary times of the day 
and year can be directly combined pixel-by-pixel. Most 
spectra were Hanning smoothed online, chiefly to combat 
the spectral rippling effects of strong narrow-band inter- 
ference somewhere in the observing band. As a result, the 
effective spectral resolution in the initial observations is 
about 150 kHz. 

The combined initial target list incorporated about 80 
sources. Despite the WSRT's impressive ability to sup- 
press external interference, by virtue of its being an inter- 
feromctric array, and having a multi-channel correlator, 
we have encountered a number of frequency intervals in 
which the level of the interference was always too high 
to obtain any useful data, even after trying to observe at 
different times, with different telescope orientations (hour 
angles), or with slightly shifted frequency bands. In the 
end, we pursued the observations for 57 sources to the 
point where the spectra showed either a line detection or 
gave a reasonable upper limit. This often involved follow- 
up observations done with narrower, offset bandwidths 
and/or twice the number of frequency channels (at the 
price of fewer correlated baselines) in the pursuit of line 
candidates. 

Bright external sources such as 3C48, 3C147, and 
3C286 were used as calibrators. Time-limited or base- 
line specific interference was first removed, before using 
them for standard complex gain and bandpass-calibration. 
No attempt was made to obtain polarisation information; 
the data correlated in the two orthogonal linear polar- 
isations were simply added after separate gain calibra- 
tion, to increase the signal-to-noise of the final spectra. 
Furthermore, since our primary interest was on line op- 
tical depths, no particular care was expended on setting 
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the overall flux density scale, but we believe that the flux 
densities shown in Figures ^ an( i are accurate to better 
than 5%. 

After external calibration, interference was removed 
from the target spectra to the extent possible. Figures 
2] and [21 show that this was not always completely suc- 
cessful: some interfering signals, particularly when their 
strength did not dominate in the observing band, could 
not be adequately discriminated by their temporal or spa- 
tial characteristics. After editing to the extent possible, 
those spectral channels without any line features were av- 
eraged, and this continuum dataset was then used in itera- 
tive cycles of self-calibration alternated with image clean- 
ing and/or model-fitting to the uv-data. The compact tar- 
gets are all unresolved to the WSRT. Characterisation of 
the continuum sources in the field was not a goal in itself, 
and was pursued only to the level needed to find adequate 
self-calibration gain factors, and to ensure that confus- 
ing components would not disturb our spectral results. 
The self-calibration complex gains found for the contin- 
uum were subsequently applied to all spectral line chan- 
nels. Final spectra were then produced by coherently inte- 
grating all phase-calibrated visibilities per spectral chan- 
nel; this method produces an optimal signal-to-noise ratio 
given that in all cases the target, at the phase centre, is 
by far the dominant source in the field. 

We find Hi absorption to be present in the spectra 
of 19 of the 57 sources (33%); these are shown in Fig.^ 
Table ^ gives the particulars of the detected absorption 
lines, determined by Gaussian profile fitting of the peak 
optical depths, r pca k, and FWHM line widths, AV; the 
column densities were computed from these assuming uni- 
form coverage of the entire radio source and a spin temper- 
ature T sp = 100 K, using N m = 1.82 x 10 18 T sp T pcak AK 
Both assumptions mean that the listed column depths are 
lower limits: the covering factor might be below unity, and 
if the H I absorption were to arise in the pc-scale vicinity 
of the AGN, conditions might be such that T sp ~ 8000 K 
(Maloney et al.[I596). For the other 38 sources, we show 
the spectra in Fig.[3 For Tabled we have used a FWHM 
of 100 km s" 1 as a plausible width for an absorption line 
in these objects (see Section l4~Tl and Fig. EJ, in order to 
be able to derive 2a upper limits to the line depths and 
H I column depths, based on the continuum flux densities 
and r.m.s. noise levels of the spectra. 

3. Detection uncertainties 

We expect that many of the optical redshifts used to cen- 
tre the initial observing bands are uncertain at the level 
of a few hundred kms" 1 ; particularly at the higher red- 
shifts, they are typically based on emission lines which can 
be quite broad, and may have a velocity centroid which is 
offset from the systemic velocity of the host galaxy. With 
a 10 MHz wide band, we have spanned a velocity range 
of approximately ±1250 km s" 1 at z = 0.2, increasing to 
approximately ±2000 km s" 1 at z = 0.85. Also note that 
an uncertainty of ±0.001 in z, which may be typical for 



the optical redshifts in the literature, corresponds to a fre- 
quency uncertainty for Hi of nearly ±lMHz at z = 0.2 
but only about ±0.4 MHz at z = 0.85. Thus, (optical) 
redshift uncertainties give the greatest (radio) frequency 
uncertainties at the lowest redshifts. In view also of in- 
terpreting the offsets discussed in Section ll~2l it would be 
valuable to obtain accurate systemic redshifts from high 
sensitivity, high resolution optical spectroscopy, for all of 
the sources we have detected in H i. Nevertheless, Fig. El 
discussed in Section l4~2l shows that in the 19 sources with 
detections, the Hi velocity distribution around the optical 
redshift tails off rapidly towards ±500 kms" 1 . We there- 
fore believe that our spanned bandwidth was adequate to 
cover the vast majority of the potential lines, and quite 
likely all of them. 



In view of the large range in flux densities amongst the 
targets, it was impractical to achieve a uniform sensitivity 
to line optical depth (opacity). It is difficult to quantify 
what this means for an individual non-detection, but the 
following limited statistical analysis is indicative. Figure|3] 
shows the individual 2a upper limits plotted along with 
the line opacities for the detected lines, while Fig.0]shows 
their cumulative distributions. Firstly, we point out that 
from Fig. [31 there appears to be no evidence for a strong 
correlation of optical depth with FWHM, which means 
that adopting a somewhat arbitrary width of 100 kms" 1 
in order to compute upper limits is reasonable (the lim- 
its would have been somewhat tighter had we adopted a 
larger width) . It is also clear that the distribution of de- 
tected line strengths is different from that of the upper 
limits: many of the limits are below most of the actual 
detections, meaning that these limits are tight enough to 
be very significant. A KS test shows that the cumulative 
distribution of the detections is different from that of the 
2a limits at the 98% confidence level. It is also clear from 
Fig. |2| however, that about half of the 2a limits leave 
something to be desired: while they rule out the presence 
of prominent absorption (roughly: above r = 1%), they 
still allow the presence of shallower lines analoguous to 
those seen in some of the weaker detections (roughly: be- 
low r = 0.5%). A significant investment in telescope time 
would be needed in most cases in order to constrain this 
potential incompleteness. 



Nevertheless, the detections that we do have are al- 
ready an important next step in studying the cold neutral 
medium in cosmologically distant radio sources. Below, we 
briefly discuss some of the line properties. In a follow-up 
paper (Pihlstrom et al. 2003), the line properties are re- 
lated to other source parameters. VLBI observations in 
order to image the neutral gas distribution are already 
being pursued for some of the most prominent absorbers 
(e.g., Vermeulen H3H1 Vermeulen et al. l2003b|l . 
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Table 1. Sources in which Hi absorption has been detected, and the line properties. z opt is the optical redshift, ^Hi.det 
the centroid frequency of the detected Hi line, Vm,dct the velocity offset of the detected line centroid from the optical 
redshift (negative means blueshifted line), the peak line optical depth, AV the line FWHM, and finally iVm the 
derived column depth. 
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GPS = Gigahertz Peaked Spectrum source, CSS = Compact Steep Spectrum source, RG = Radio Galaxy, 
BLRG = Broad Line Radio Galaxy. 

The radio band centre frequency used during the observations sometimes differed from that predicted by the optical 
redshift, for example in order to avoid RFI or to pursue tentative lines. 
Derived from 7V H i = 1.82 x 10 18 T S pT pcak AV cm" 2 , using T sp = 100 K. 

The detection of Hi absorption in 3C 216, an integral part of our survey, was first published in Pihlstrom et al. (1999 ) 
Redshift based on Ca H+K and G-band stellar absorption features in an unpublished spectrum obtained by some of us 
(RCV and GBT) at the Palomar Observatory 200" Hale Telescope. 



4. Line Properties 

4.1. Line widths 

Figure shows the distribution of the FWHM of the 
Gaussians fitted to the absorption lines. This plot con- 
firms that H I lines in these sources often are broad, with 
a FWHM of 100-200 km s _1 or even more. We find that 
the primary lines (i.e. the deepest Gaussians, with the 
highest optical depth in each source) are probably some- 
what broader than the secondary lines; the mean FWHMs, 
of 156kms _1 and 72kms~ 1 , respectively, differ at a confi- 
dence level of 97%, according to a Student's i-test. 



Often, Hi absorption is interpreted as being due to 
a circumnuclear disk; for example VLA 21 cm Hi ab- 
sorption imaging of a sample of Seyfert and starbursts 
displayed almost exclusively sub kpc scale rotating disks 
aligned with the outer galaxy disk (Gallimore et al. I1999|) . 
Data for radio loud sources are also consistent with disks 
on sub kpc scales, for example in Cygnus A fGonwav ll999[) 
and 1946+708 (Peck & Tavlor l2T)0T)> . Those examples all 
have line widths of the order of 100 km s -1 , comparable 
to the detections reported in this paper. If the wide lines 
arise on sub kpc scales close to the nucleus, the narrower 
lines could be additional absorption occurring on larger 
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Fig. 1. Spectra of sources with detected Hi absorption. The fitted Gaussian profiles are also shown, and features due 
to radio frequency interference are marked 'RFF on those plots. Velocities are with respect to the optical redshifts. The 
grey horizontal bands represent the la standard deviation from the fitted continuum level. The grey vertical regions 
indicate the outer 25% of the total bandwidth. 



scales within the galaxy, provided that the line of sight 4.2. Inflowing and outflowing gas ? 

does not cross a substantial amount of differential rota- 
tion. Alternatively, the narrower lines could be from sin- It was pointed out by van Gorkom et al. (1989) that in all 
gle cohesive clouds, with only moderate internal velocity 4 out of 29 nearby radio galaxies (z < 0.13) with detected 
dispersion, as is observed in 3C236 (Conway & Schilizzi Hi absorption, and also in 4 other nearby systems, the 

Hi velocity is positive, i.e. infalling with respect to the 
systemic velocity of the galaxy. The radio structure (or, 
in the case of 3C236, at least the absorbed part) of those 
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Table 2. Sources in which Hi absorption has not been detected, with upper limits to the lines. z opt is the optical 
redshift, t 2(J the limit to the optical depth for a line of width AV = 100 km s -1 centred at the optical rcdshift, and 
finally N H \_ 2a the corresponding limit on the column depth. 
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1 GPS = Gigahertz Peaked Spectrum source, CSS = Compact Steep Spectrum source, RG = Radio Galaxy, 
RQ = Radio Quasar, BLRG = Broad Line Radio Galaxy, FRII = Fanaroff-Riley Type 2, 

EORG = End On Radio Galaxy and EORQ = End On Radio Quasar. 

2 The radio band centre frequency used during the observations sometimes differed from that predicted by the optical 
redshift, for example in order to avoid RFI or to pursue tentative lines. 

3 Derived from JV H i, 2ct = 1.82 x 10 18 T sp r 2CT AV cm~ 2 , using T Bp = 100K and AV = 100 km s -1 . 

4 Redshift based on Mgll, [Oil], and [OIII] emission lines in an unpublished spectrum obtained by some of us (RCV 
and GBT) at the Palomar Observatory 200" Hale Telescope. 
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Fig. 2. Spectra of sources without detected Hi absorption. The horizontal lines shown represent the ler standard 
deviation from the fitted continuum level. Features due to radio frequency interference are marked 'RFF. Velocities 
are with respect to the optical redshifts, positive corresponding to infall. The grey vertical regions indicate the outer 
25% of the total bandwidth. 
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Fig. 2. (cont). Spectra of sources without detected Hi absorption. The horizontal lines shown represent the la 
standard deviation from the fitted continuum level. Features due to radio frequency interference are marked 'RFI'. 
Velocities are with respect to the optical redshifts, positive corresponding to infall. The grey vertical regions indicate 
the outer 25% of the total bandwidth. 



nearby galaxies with Hi absorption is very compact, and 
many of these objects are now classified as GPS or CSS 
sources. Thus, they potentially match well with the types 
of sources that constitute the majority of our sample. 

As Figure shows, we now find Hi lines at both pos- 
itive and negative velocities with respect to the optical 
redshift. Also, while typical velocities in the lower red- 
shift sources of van Gorkom (1989) were in the range to 
+200 kms -1 , with only 1 value around +400 kms -1 , 10 of 
our 19 main (highest opacity) lines have a velocity more 
than 200 kms -1 different from the optical redshift; the 
main line velocities range between v = — 1420 kms -1 and 
v = +318 kms -1 . While, as already discussed in Section[SJ 
inaccuracies in the optical redshift may play a role in caus- 
ing some of the scatter, we believe the larger H I velocities 



are likely to be significant. We see an intriguing hint that 
the velocity distribution may be skewed towards outflow 
rather than infall. This effect can be stated more dra- 
matically than it appears through the binning of Fig. |SJ 
13 of the main line velocities are negative as compared 
to 6 positive; of these, 7 have v < —200 kms -1 as com- 
pared to 3 having v > 200 kms -1 , and of these, 3 have 
v < —600 kms -1 as compared to none at v > 600 kms -1 . 
The mean centroid velocities are v = —223 kms -1 for the 
main lines and v — — 58kms _1 for the secondary lines, re- 
spectively. However, Student's t tests at a confidence level 
of 95% show that the mean velocities of the main lines and 
the secondary lines indeed do not differ significantly from 
zero, and also not from each other, possibly only because 
the sample size is too small. 
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Fig. 3. Distribution of line strengths and widths. The 
2(7 non-detection limits are shown as vertical marks at 
100 km s™ 1 ; any Hi absorption in these objects probably 
has an opacity smaller than (to the left of) the mark. 
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Fig. 4. The cumulative distribution of the opacity of the 
detections (solid line) and the non-detections (dashed 
line), as discussed in Section 
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Fig. 5. Number distribution of line widths, in 40kms~ 1 
wide bins, as discussed in Section 14.11 The main lines 
(largest optical depth in each source) and the secondary 
lines are indicated separately. 
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Fig. 6. H I line velocities compared to the optical redshifts, 
in bins of width 200 kms -1 . The main lines (largest optical 
depth in each source) and the secondary lines are indicated 
separately; their distributions are discussed in Section f4.2l 



Thus, we find the opposite of the results obtained for 
nearby radio sources by van Gorkom et al. ( 1989) in both 
the sign and the magnitude of the velocities. Speculatively, 
we think that in these higher luminosity sources, inter- 
action between the radio jets and the surrounding inner 
galactic medium could lead to significant motions of the 
gas, such as found in 3C236 by Conway & Schilizzi (2000). 
Other evidence for such interactions is often seen as the 
so-called radio-optical alignment effect. The kinematics of 
the aligned gas seen in CSS sources shows velocity offsets 
of a few hundred kms -1 (both positive and negative) with 
respect to the nuclear velocity (e.g., O'Dea et al. [5002 ) . 
Alternatively, in some cases we may be looking at the 
line-of-sight rotational motion of a neutral component of 
broad-line gas, seen in front of a parsec-scale radio source, 
as has been discussed for lower redshift galaxies for exam- 
ple by Morganti et al. 



5. Summary 

Hi absorption searches have been carried out with the 
WSRT for a sample of mostly compact radio sources with 
redshifts out to nearly 2 = 1, giving observing frequen- 
cies which fall in the UHF-high band (725-1200 MHz). 
Out of 57 sources for which adequate spectra could be ob- 
tained, 19 are found to have associated H I in absorption, 
with opacities ranging from 16% down to 0.2%; below 1% 
opacity, our statistics may well be incomplete. The typi- 
cal line width is ~ 150 kms -1 , but the lines can also be 
several times narrower or wider. Column densities ranging 
from 1 x 10 19 to 3.3 x 10 21 are obtained by assuming uni- 
form coverage of the radio source and a spin temperature 
T sp = 100 K. In contrast to nearby radio sources, we find 
not only positive but even more negative Hi velocities, 
up to more than 1000 kms -1 (as compared to the optical 
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redshifts). Perhaps, inflow and outflow can both occur, or 
there is interaction with the radio source, or we are seeing 
rotating broad-line gas. 
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